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Abstract 

Based on a faithful representation of the heavy tail multivariate distribution of asset returns 
introduced previously (Sornette et al., 1998, 1999) that we extend to the case of asymmetric 
return distributions, we generalize the return-risk efficient frontier concept to incorporate the 
dimensions of large risks embedded in the tail of the asset distributions. We demonstrate that 
it is often possible to increase the portfolio return while decreasing the large risks as quantified 
by the fourth and higher order cumulants. Exact theoretical formulas are validated by empirical 
tests. 
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1 Introduction 



One of the most fundamental tenet of eeonomie theory and practice is that returns above the 
so-caUed riskless rate come with increased risks. This is the basis of Markovitz's portfoho theory 
(e.g. Markovitz, 1959) and of the CAPM (e.g. Merton, 1990). Reciprocally, investors want to 
be compensated for taking risk, that is, they want to earn a return high enough to make them 
comfortable with the level of risk they are assuming. It is thus a fundamental premise of efficient 
markets that "one cannot have both the cake and eat it too", i.e. one cannot increase the return 
and lower the risk at the same time. This result stems simply from the linear (resp. quadratic) 
dependence of the average return (respectively variance) of a portfolio return on the weights of its 
constituting assets leading to a parabolic efficient frontier in the return-risk diagram. 

In the real world, the variance of portfolio returns provide only a limited quantification of incurred 
risks, as the distributions of returns have "fat tails" (e.g. Lux, 1996, Gopikrishnan et al., 1998, Lux 
and Marchesi, 1999) and the dependences between assets are only imperfectly accounted for by the 
correlation matrix (e.g. Litterman and Winkelmann, 1998). Value-at-Risk (e.g. Jorion, 1997) and 
other measures of risks (e.g. Artzner et al., 1996, Sornette, 1998, Bouchaud et al., 1998, Sornette 
et al., 1998, 1999) have been developed to account for the larger moves allowed by non-Gaussian 
distributions. 

Here, we generalize our previously introduced representation of the heavy tail multivariate distri- 
bution of asset returns (Sornette et al., 1998, 1999) to the case of asymmetric return distributions. 
We calculate theoretically and test empirically the cumulants of a portfolio and generalize the 
return-risk efficient frontier concept to incorporate the dimensions of large risks embedded in the 
tail of the asset distributions. We demonstrate the novel remarkable result that it is often possible 
to improve on the optimal mean-variance portfolio by increasing the return while decreasing the 
large risks quantified by the fourth and higher order cumulants. This is related to and generalizes 
our previous rigorous result (Sornette et al., 1998, 1999) that minimizing the variance, i.e. the 
relatively "small" risks, often increases larger risks as measured by higher normalized cumulants 
and the Value-at-risk. Thus, putting the emphasis on the risk quantified by the volatility can be 
both misleading because large risks are still looming and in addition damage profitability. 



2 The asymmetric modified Weibull distribution 

In order to make our approach concrete, we assume that price returns Sx are distributed according 
to the following probability distribution function (pdf) 

P{6x) = -^^^ \dx\'-T-^ exp(^-(^^y^^ forO<fe, (1) 

= ^^M--' -p(-(f^)") P) 

Q is the probability for observing a positive return, the x's are the characteristic returns and the 
exponent 7's control the fatness of the pdf tails, which can be different for positive and negative 
returns. 
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For Q = 1/2, x+ = X- aiid 7+ = 7_, we recover the symmetric modified Weibull pdf studied by 
Sornette et al. (1998, 1999) and the special case 7+ = 7_ = 2 recovers the standard normal law. 
The case when the exponents 7 are smaller than one corresponds to a "stretched" exponential with a 
tail fatter than an exponential and thus much fatter than a Gaussian, but still thinner than a power 
law. Stretched exponential pdf 's have been found to provide a parsimonious and accurate fit to the 
full range of currency price variations at daily intermediate time scales (Laherrere and Sornette, 
1998). This stretched exponential model is also validated theoretically by the recent demonstration 
that the tail of pdf 's of products of a finite number of random variables is generically a stretched 
exponential (Frisch and Sornette, 1997), in which the exponent 7 is proportional to the inverse of 
the number of generations (or products) in a multiplicative process. 



3 Nonlinear change of variable 

Let us pose 



y+ = ((5x)^+/2 for 6x > , (3) 
y_ = -|<5xp-/2 for 6x < . (4) 



Inversely, we have 
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6x = y'^ for 6x > , with g+ = — , (5) 

2 

5x = — ly-l*^" for 5x < , with g_ = — . (6) 

7- 

The change of variable (^,^ from 6x to y leads to a Gaussian pdf for the y-variable defined in each 
semi-infinite domain: 
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P(y_) = ^^^^ exp ( ) , where al = ^xl" ■ (8) 



Using a maximization entropy principle, one can then show (Sornette et al., 1998, 1999) that the 
correlations between the y variables of different assets provide the most efficient and parsimonious 
multivariable representation. This transformation has also been used for the analysis of parti- 
cle physics experiments (Karlen, 1998) and much earlier for the treatment of bivariate gamma 
distributions (Moran, 1969). It can also be viewed as a concrete implementation of the copula 
representation of dependence between assets (e.g. Embrechts et al., 1998, 1999). Generalizations 
to other non-Gaussian pdf's are discussed in Sornette et al. (1999). 

We have made empirical tests on three assets, using annualized daily returns of stock prices of 
Chevron (CHV) and Exxon (XON) in the period Jan. 1970 - Mar. 1999, and of the Malaysian 
Ringit (MYR) against the US dollar in the period Jan. 1971 - Oct. 1998. The CHV-XON pair 
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is among the most strongly connected group of stocks in the S&P 500 index while the Malaysian 
Ringit is essentially uncorrelated to the Chevron and Exxon stocks. These extreme cases allow 
us to test the influence of correlations. Especially for strongly correlated stocks, we have shown 
(Sornette et al., 1998) that a change of variable like Eq. leads to a covariance matrix which is 
much more stable compared to the usual covariance matrix. 

Fig. |l| shows in a log-log plot the y(r) transformation (^j^ calculated from the empirical positive 
and negative returns of the Chevron and Exxon stocks and for the Malaysian Ringgit against the 
US dollar (MYR) . Assuming that price returns are distributed according to an asymmetric modified 
Weibull ([l|,§), the slope of the y(r)-plot gives for large |r|-values the exponents 7+/2 and 7-/2. The 
positive and negative returns of each asset are seen to have almost the same slope for large r values, 
and consequently we will assume for each asset that 7+ = 7_ = 7 in the sequel. The linearity of 
the y(r) plots for large r values show that the large tails of the pdf's are indeed to a very good 
approximation distributed according to a modified Weibull distribution Eq. (^,0), with 7 ~ 1.4 
(CHV), 7 ~ 1.2 (XON) and 7 « 0.62 (MYR). For small and intermediate r values, the y{r) curves 
have a slope close to 1 (indicated by the y = r line), which means that small and intermediate 
returns are distributed according to a Gaussian distribution. Because of the finite resolution of the 
data (the data has a lower bound for the return), y{r) approaches a constant value for the smallest 
values of r. 



4 Portfolio theory for the diagonal case 

In this short letter, we present the theory for the diagonal case where assets are uncorrelated. 
This is already sufficient to illustrate the most important results. Especially in the case of fat tails 
(exponents c < 1), correlations are less important than a precise determination of the tails (Sornette 
et al., 1998). We will however present some empirical tests with uncorrelated and with correlated 
assets, in order to illustrate the importance of correlations. Sornette et al. (1999) treat the case 
of correlated assets with symmetric distributions with the same exponent 7. Generalization to the 
asymmetric case and with different exponents 7 will be reported elsewhere. 

The discrete time estimation of the returns 6xi{t) are 6xi{t) = 6pi{t)/pi{t) = {pi{t + l) — pi{t)) /pi{t) , 
where Pi{t) is the price of asset i at time t. The total variation of the value of the portfolio made 
of N assets between time t — 1 and t reads 

N N 

6S{t) = ^ Widpiit) = J2wi6xi{t) , (9) 

i=l 1=1 

where Wi is the number of shares invested in asset i and Wi = WiPi is the weight in capital invested 
in the ith asset at time t in the portfolio. We will assume normalization, i.e. J2iLi Wi = 1, thus 
leading to a dynamical reallocation of the assets in the portfolio. 

The expression (|9|) can be expressed in terms of the variables yiS defined by (^,^) as follows 

N 

<55(t) =^u;, e.lyej"'' , (10) 

i=l 
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where ej is the sign of 6xi. All the properties of the portfolio are contained in the probability 
distribution Ps{SS{t)) of 6S{t). We would thus like to characterize it, knowing the multivariate 
distribution of the 6xi's (or equivalently the multivariate Gaussian distribution of the y^'s) for the 
different assets. The general formal solution reads 



Psi6S) =cf[([ dy,) 6-^2 y'v-'y sfdSit) -j^w^e^ly^^'^A . 

i=i V"' / V i=i / 

Taking the Fourier transform Ps{k) = dSS Ps{SS) e''^^^ of ([l|) gives 

Ps{k) = {[([ dy,) e~-2 y'^-V+^fc Eti ^d^/sl"' . 
i=i ^ 

Using the explicit expression of the form of the distributions (0,^ , we get 



Ps{k) = n 

i=l 
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Expanding the exponential exp {ikwi\yi\'^'^) in powers of its argument, we get 

{ikwi) 
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and r is the Gamma function. Replacing in (14), we obtain 
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Psik) = n 



i=l 
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where 



(11) 



(12) 
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(14) 



(15) 



(16) 



(17) 

For symmetric distributions with qi+ = qi^, i.e. 7^+ = 7i_, C7j+ = cjj_ and Qi = 1/2, we retrieve 
our previous result (Sornette et al., 1999) that all the odd order terms in the sum over m cancel 
out : 



N 

Psik) = n 

i=l 



.n=0 



(2n)!^ 



(18) 
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The expression X/m^o ^^^ni\ — ^^ii"^) ( |16D is similar to the expansion of a characteristic function 
in terms of moments. We need to get the corresponding expansion in terms of cumulants, i.e. find 
the coefficients Cn such that 

£ M.(n.) = exp rf M! . (19) 

m=0 \?i=l / 

By identifying the same powers of k term by term, we get the cumulants. Then, using the product 
in (16) of the exponentials from i = 1 to A^, we obtain the cumulants of the portfolio distribution 
as 



N 



Cl 



Y^WiMiil) , (20) 

i=l 
N 

J2wf [M,i2) - M.ilf) , (21) 
1=1 

N 

Y,wf{Mii3)-3M,il)Mi{2) + 2Miilf) , (22) 

i=l 
N 

C4 = (Mi(4)-3M,(2)2-4Mi(l)Mi(3) + 12Mi(l)2Mi(2)-6Mi(l)4) , (23) 

1=1 

N , 

C5 = E """l ^^(5) - 5Mi(4)Mi(l) - 10Mi(3)M,(2) + 20Mi(3)Mi(l)2 + ?>^Mi{2f MiiX) 



C2 



C3 



i=l 



60Mi(2)Mi(l)3 + 24Mi(l)M , (24) 



N 



C6 = E"^n ^^^(6) - 6Mi(5)Mi(l) - 15Mi(4)Mi(2) +30Mi(4)Mi(l)2 - X^MiiZf 
i=i ^ 

+120Mi(3)Mi(l)Mi(l) - 120Mi(3)Mi(l)3 + 30Mi(2)3 

-21^Mi{2fMi{\f + 360Mi(2)Mi(l)^ - 120Mi(l)M . (25) 



Higher order cumulants are obtained by using the formulas given for instance by Stuart and Ord 
(1994). The first cumulant c\ provides the average gain l^bS) and the second cumulant ci is the 
variance of the portfolio gain. The higher order cumulants as well as the excess kurtosis k = C4/C2 
quantify larger risks occurring with smaller probabilities but larger impact. 

Fig. § presents a comparison of the empirical determined c„'s and those determined from the 
equations ( p0[]2^ ), for a portfolio constituted of the Malaysian Ringgit (MYR) and the Chevron 
stock (CHV). This choice is made because MYR is essential uncorrelated to CHV and the above 
calculation should thus apply directly. For an extension of the theory to correlated assets, see 
Sornette et al. (1999). To perform the empirical test shown in figure |2|, we first determined the 
exponents 7+ = 7_ = 7 from a regression of the linear parts of the y{r) functions for large values 
of |r| shown in figure 1. We then use these 7's to estimate the coefficients Xi+^Xi- from the 
empirical averages Xi± = {{^x±)"''-)±. The notation ( )± represents an average taken with respect 
to positive/negative returns of the data. The asymmetric weight parameter Qi is determined from 
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the asset i as the ratio of the number of positive returns over the total number of returns. The error 
bars shown in the figure are determined from the observation that the main source of error comes 
from a mispecification of the tail exponent 7's and we assume conservatively an error of ±0.05 
on the 7 values. Fig. |2| shows a very good agreement between theory and the direct empirical 
determination of the cumulants. There is some discrepancy for the third order cumulant C3, which 
reflects our simplification to use symmetric tails with 7+ = 7_ = 7 in our calculations Eq. 
As a consequence, the sole contribution to the odd-order cummulants stems from the difference 
between Xi+ ^^id Xi- ^-^d between Qi and 1/2. An additional asymmetry in the shape of the 
tail captured by 7+ 7^ 7-, however small, can easily make the agreement adequate between the 
theoretical and empirical C3. We have chosen not to incorporate this additional complexity in order 
to keep the number of degrees of freedom as small as possible. The even-order cumulants and the 
excess kurtosis k are much less sensitive to the asymmetry in the exponents 7_|_ , 7_ . 

The portfolio with minimum variance C2 has the optimal weight wi = 9.5%, where the index 1 
stands for the Chevron stock, i.e. the weight W2 = 1 — wi of the Malaysian Ringgit is 90.5%. In 
comparison, the portfolio with minimum fourth cumulant has an investment ratio of wi = 38% in 
Chevron and W2 = 62% in the Malaysian Ringgit. It is clear that the minimum variance portfolio 
has a rather large fourth cumulant, i.e. minimizing the small risks quantified by the second order 
cumulant comes at the cost of imcreasing the largest risks quantified by the fourth order cumulant 
(Sornette et al., 1998, 1999). 

Fig. ^ illustrates another even more interesting phenomenon. We compare the daily returns and 
the cumulative wealth of two portfolios. The first ci — C2 portfolio has a minimum variance C2 
(Chevron weight wi = 0.095 and Malaysian Ringgit weight W2 = 0.905). The second ci — C4 
portfolio has a minimum fourth-order cumulant (Chevron weight wi = 0.38 and Malaysian Ringgit 
weight W2 = 0.62). The horizontal dotted lines in the daily return plots are the maximum values 
sampled for the returns of the ci — C4 portfolio. Notice that the daily returns of the minimum 
variance portfolio exceeds these bounds. This illustrates vividly that, while most of the time the 
fluctuation of the returns are smaller for the ci — C2 portfolio, fluctuations with larger amplitudes 
and thus larger risks are observed in this minimum variance portfolio: again, minimizing small 
risks can lead to a dangerous increase of large risks (Sornette et al., 1998, 1999). Furthermore, the 
cumulative wealth of the ci — C2 portfolio with wi = 0.095 is drastically inferior to that accrued in 
the ci — C4 portfolio with wi = 0.38. In other words, you can have your cake and eat it too: decrease 
the large risk (those that count for the safety of investment houses and for regulatory agencies) and 
increase the profit! This example illustrates how misleading can be the focus on the variance as a 
suitable measure of risks and how limited is the use of standard portfolio optimization techniques. 
Not only they do not provide a suitable quantification of the really dangerous market moves, in 
addition they miss important profit opportunities. 

Fig. ^ is the same as Fig. ^ for a portfolio constituted of the Exxon and the Chevron stocks. Due to 
the very large correlation between the two assets, the departure between theory and experiments 
is a measure of the importance of correlations that have been neglected in the above formulas, 
expecially in this case where the exponents 7 for the pdfs of the two stocks are relatively large 
around 1.4 and 1.2 respectively, i.e. the pdf tails are relatively "thin". This constitutes a worst- 
case scenerio for the application of the above theory that is best justified for exponents 7 < 1 (recall 
that the standard Gaussian regime corresponds to 7 = 2). Nothwithstanding this limitation, the 
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results conform qualitatively to our previous discussion: the best variance gives a substantially 
larger risk for large moves and the return is sub-optimal. 



5 Efficient Portfolio Frontiers 



Based on our previous calculation, it is straightforward to construct the optimal mean-variance 
portfolios from the knowledge of the cumulants ci and C2 as a function of the asset weights Wi. 
Similarly, we introduce the optimal ci — C4 portfolios. 

For a given mean return ci, the portfolios that minimize the risks expressed through C2 given by 



(26) 
(27) 



Eq. (21) or by C4 given by Eq. (|23|))) are determined from the conditions 

d 



d 



C2 - AlCi - A2^a;j 

i 

C4 - AiCi - A2 ^ 



, 
, 



where the ujj denote the weights for an optimal portfolio. From the normalization condition 

E^i = l' (28) 

i 

one of the Lagrange multipliers among Ai, A2 can be eliminated. Let us define cn such that the 
expressions (20,^,^) read 



i 
i 

C4 = Eu;fc4i . 

i 

The efficient frontier for the mean-variance ci — C2 porfolios is given by: 
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-}-{A-B'/D) + B/D 
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Varying Ai then traces out the efficient frontier. Likewise the efficient frontier for the ci — C4 
portfoUos is given by: 



Cl = ^c^*cl, , (38) 

i 

C4 = ^(w*)^c4i , with (39) 

i 

^* ^ I , | (clj - A2) ,4/3 /.nN 

' E.±l(cl^-A2)/(4c4,)|i/3 I 4c4, ' ' ^ ^ 

with + if clj > A2 and — otherwise. 

Fig. ^ shows the efficient frontiers for portfohos constituted of the three assets CHV-XON-MYR. 
The hnes are derived from the theoretical prediction given by Eq. (|27| ) using the exponents de- 
termined from Fig. |l|. The sohd hne shows the mean-variance efficient frontier normahzed to the 
minimum variance and the dotted hne shows the ci — C4 efficient frontier normahzed to the minimum 
fourth-order cumulant determined from the theory assuming no correlations between the assets. 
The -|- (resp. o) are the empirical mean-variance (resp. ci — C4) portfolios constructed by scanning 
the weights wi (Chevron), W2 (Exxon) and W3 (Malaysian Ringgit) in the interval [0,1] by steps 



of 0.02 with the condition of normalization (28). Both family define a set of accessible portfolios 
and the frontier of each domain define the corresponding empirical efficient frontiers. Note that by 
allowing negative weights (short position), the domains within the parabola are progressively filled 
up, corresponding to accessible portfolios with "short" positions. 

The agreement is not good quantitatively between theory and empirical tests due to the strong 
correlations between Chevron and Exxon which is neglected in the theory (see figure (§)). How- 
ever, there is good qualitative agreement: the theory and empirical tests are essentially translated 
vertically, with the same characteristics. The most important feature is that the ci — C4 portofolio 
with minimum fourth-order cumulant (small "large risks" ) has a significantly larger return ci than 
the portfolio with the minimum variable. For instance in the historical data, the return for the 
minimum variance occurs for wi = 0.032, W2 = 0.084, W3 = 0.884 for which the mean annualized 
return is ci = 3.1% and the fourth-order cumulant is C/^/c/^min = 2.22, i.e. more than twice the 
minimum possible value. The minimum of C4 is reached for wi = 0.292, 11)2 = 0.084, 1^3 = 0.624 for 
which the mean annualized return is ci = 7.2%, i.e. more than double the return for optimal the 
mean-variance portfolio. Its variance is C2/c2min — 1.T3 which is a relatively moderate increase of 
"small risks" . The results presented here can be easily generalized to higher cumulants with similar 
conclusions. 
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Figure 1: ?/(r)-transformation defined by equations for the period from January 1971 to oct. 
1998. + corresponds to positive returns and o to negative returns. The daily returns r are expressed 
in annuahzed percentage, a) Chevron stock (CHV), b) Exxon stock (XON), c) Malaysian Ringgit 
against US dollar (MYR). 
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Figure 2: Comparison of the empirically determined cumulants c„ and excess kurtosis k (fat solid 
line) to the theory Eq. (|20|-p3|) (thin solid line) using the exponents determine from Fig. || for a 
portfolio constituted of the Malaysian Ringgit and the Chevron stock. The cumulants are plotted 
as a function of the asset weight wi, where the index 1 corresponds to CHV, with the normalization 
vji + W2 = 1. Thus, the weight of the Malaysian Ringgit is t(;2 = 1 — if i • The error bars shown are 
obtained assuming an uncertainty in the determination of the exponents 7^ = 7^ it 0.05. 
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Figure 3: Annualized daily returns (in percent) and cumulative wealth (starting with a unit wealth 
at time zero) for the two portfolios corresponding to the minimum variance with Chevron weight 
wi = 0.095 and minimum fourth-order cumulant C4 with Chevron weight u^i = 0.38, determined 
from figure |2[ 
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Figure 4: Same as figure ^ for a portfolio constituted of the Exxon and the Chevron stocks. The 
cumulants are plotted as a function of the Chevron weight wi and the weight of the Exxon stock 
is 102 = 1 — Wl. 
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Figure 5: Efficient frontiers for the three-asset portfolio CHV-XON-MYR derived from theory 
Eq. (^) using the exponents 7j's determined from Fig. (|l]). The solid line shows the mean- variance 
efficient frontier normalized to the minimum variance and the dotted line shows the ci — C4 efficient 
frontier normalized to the minimum fourth-order cumulant determined from the theory assuming 
no correlations between the assets. The -|- (resp. o) are the empirical mean-variance (ci — C4) 
(resp. ci — C4) portfolios constructed by scanning the weights wi (Chevron), W2 (Exxon) and 
(Malaysian Ringgit) in the interval [0, 1] by steps of 0.02 while still implementing the condition of 
normalization ( p8|) . Both family define a set of accessible portfolios excluding any "short" positions 
and the frontier of each domain define the corresponding empirical frontiers. 
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